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Abstract 


A direct derivation is given of a formula for the normalized asymptotic variance 
parameters of the boundary local times of reflected Brownian motion (with drift) on 
à compact interval. This formula was previously obtained by Berger and Whitt using 
an M/M/1/C queue approximation to the reflected Brownian motion. The bivariate 
Laplace transform of the hitting time of a level and the boundary local time up to that 
hitting time, for a one-dimensional reflected Brownian motion with drift, is obtained 
as part of the derivation. 
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In [2], Berger and Whitt studied diffusion approximations to a rate-control throttle 
and the G/G/1/C queue. These diffusions are one-dimensional reflected Brownian 
motions (with drift) on the interval [0, C], where C < œ. Amongst the formulas derived 
in [2] were expressions for the normalized asymptotic variance parameters of the 
boundary local times (or boundary regulator processes) of these reflected Brownian 
motions (see (4.14) of [2] and (1)-(6) below). These normalized asymptotic variance 
parameters are of interest because they are the squared coefficients of variation (SCVs, 
variance divided by the square of the mean) for the constituent 1.i.d. random variables in 
a renewal process approximation to the overflow processes in the rate-control throttle 
and queueing models considered by Berger and Whitt [2]. Berger and Whitt obtained the 
expressions for these asymptotic variance parameters by first showing that they are 
limits of quantities for M@/M/1/C queues and then explicitly computing the M/M/1/C 
quantities ({2], Theorem 3.1) and their limits ([2], (4.36) of Theorem 4.5). Ward Whitt 
asked the author whether one could not perhaps calculate the expressions for the 
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parameters in (4.14) of [2] directly from an analysis of the reflected Brownian motion. In 
this note it is shown that this is indeed possible. The objective is to compute a quantity 
(T defined in (5) below) which is made up of moments and joint moments of the hitting 
time of a level and the boundary local time up to that hitting time for a one-dimensional 
reflected Brownian motion with drift. The moments and joint moments are computed 
by differentiation of a bivariate Laplace transform. A formula for the latter is obtained 
by a deft choice of a suitable martingale. The symbolic differentiation of the Laplace 
transform was performed using the computer software package Mathematica on a SUN 
3/60. 

We begin by defining the reflected Brownian motion (RBM) under consideration. Fix 
C>0. Let X be a one-dimensional Brownian motion with drift „ER, variance 
parameter a” > 0, and suppose X(0) = 0. Define A = o { X(s):0 <s St} forall t Z0. A 
process will be called adapted if it is adapted to {A}. There is a unique pair of 
non-decreasing, continuous, adapted, one-dimensional processes (L, U) such that 

(i) Z(t) = X(t) + L(t) — U(t) E[0, C] for all t 2 0, 

(i1) L can increase only when Z is at 0, 

(iii) U can increase only when Z is at C. 

Existence and uniqueness follow easily from the more general work of Lions and 
Sznitman [6], or a ‘bare-hands’ proof can be given by successive application of the one- 
sided reflection mapping ([3], Chapter 8) using a sequence of stopping times. The process 
Z is referred to as a (£, o°) reflected Brownian motion on [0, C] with starting point at the 
origin. The processes L and U are called the local times of Z at the origin and C 
respectively. Now, define t = inf{t 20: Z(t) = C}, t’ =inf{t 271: Z(t) =0} —t. As 
claimed in ({2], §4.2), by a regenerative analysis and the central limit theorem for 
renewal processes, one has 


(1) a=lim EO — EC + 
t= fÍ E[t + T] 


. U(t) £E[U(t +7’)] 
m-_ = MlM 


(2) b =li 
t=o t Eft +1] 
Var(L(t Var(U 
(3) lim ae =o7 and lim _—_ = Gj, 
[+2 t— œ 


where oż , oł, are called the asymptotic variance parameters of L and U respectively, and 
the associated normalized asymptotic variance parameters I; , I, defined by 


(4) T; =a 'a? and Ty =b-'o}, 


are the squared coefficients of variation for the i.i.d. random variables constituting 
renewal process approximations to the overflow processes in the models considered by 
Berger and Whitt ((2], (4.13)). In fact, it turns out that I; = Ty (see [2], (4.14) or Remark 
2 following Theorem 2 below), so we shall henceforth only consider I, and denote it 
simply by I. By a regenerative analysis ({2], (4.45)) and the fact that L does not increase 
on [t,t + 7’), 
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t+ t’)E[L(t)]\2 
The main object of this note is to verify formula (4.14) of [2], namely, 
2C7/3, ifu = 0 
(6) "= }2(1 —exp(26C))+40C exp(6C) 
MC if u #0, 


— 0(1 — exp(0C))? ? 


where 0 = 2u/0°. In order to do this, the first and second moments of t, t’ and L(t) need 
to be computed, as well as the first moment of tL (rt). 

Note that if o > 0 and Z is a (u/c, 1) reflected Brownian motion on [0, C/a], then by 
Brownian scaling and the representation of the form (i) for Z, we have Z =aZisa (u, o?) 
reflected Brownian motion on [0, C] and the local time L of Z at the origin is equal to oÊ 
where L is the local time of Z at the origin. It follows that it suffices to verify (6) for 
a? = 1. Thus we henceforth assume o? = 1. 

The following theorem is the key to our computation of I. The constituent moments 
required to compute T, that can be derived from (7), are displayed in Theorem 2, 
towards the end of this paper. These moments also enable one to compute all of the 
asymptotic quantities defined in (1)-(4) above. 


Theorem 1. For each a, B 2 0 we have 


] 
if u = 0,8 =0 
l +a&C me 
] 
E ifu =0, B #0 
(T) Elexp( — aL(t) — Bt)] = 4 cosh(yC) + i sinh(yC) 
uC 
7 = ifu ¥ 0, 
cosh(yC) + S sinh(yC) 
y 


where y = (u? + 2B)". 
Proof. Consider the function f: R X R? —R defined by 


K(x, 1, t) 
(8) 
= exp( — u(x + /))[y cosh(y(x + D) + (a + w)sinh(y(x + /))Jexp( — al — Bt). 


It can be readily verified that f satisfies: 
(9) bfx tf +h =0 


and 
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(10) f(x, -—x,t)=90 for all (x, ER X R,. 
Then by It6’s formula we have for B, =X, — us, 
Xxx | er [A T) — f(0, 0, 0) 


(11) = [T RA, Le, 94B, + | 7 AA, Les s)dL, 
0 0 
+ [T OSa HUR + INK, La, Dds. 
0 


Since U does not increase prior to the time t and L(- AT) can increase only when 
Z(-At)=X(- at) +L(- AT) 1s at the origin, it follows that the second integrand in (11) 
can be replaced by f,(X,, — X,, s), and then using (9)-(10) we see that 


(12) A\Xines Lins LAT) = f(O, 0, 0+ | LUG, Le, 9)4B, 
0 


which is a local martingale. In fact, since Z, = X, + L, E[0, C] for OSs tT and 
aL, + Bs 2 0, it follows from (8) that the left member of (12) is bounded by a constant 
independent of t, and hence that it is actually a martingale (cf. Proposition 1.8 of [3]). 
Indeed, since tT < œ a.s., after taking expectations in (12) and letting t — œ we obtain by 
bounded convergence that E[ f(X., L,, t)] = f(0, 0, 0). That is, 


E[e~““(y cosh(yC) + (a + w)sinh(yC))e~ 2-4] =y. 


Note that the argument in the expectation above is always positive for y # 0 since 
(|u|/y)tanh(yC) < 1 in this case. Thus for y # 0, formula (7) is verified. If y = 0, then 
u = 0, $ = 0, and the formula in the first line of (7) can be obtained by letting y (hence £) 
tend to 0 in the second line of (7). 


Remark i. For the case a> 0, uw = 0, formula (7) can alternatively be obtained from 
a formula given in ([7], Exercise (VI.4.9), part 4°) for u =0, by using a Girsanov 
transformation and the fact that (Z,L) = (|X|, L°) when u =0, where Z(t) = 
X(+ L(t), L(t) = — miny<,<, X(s), and L°? is the local time at the origin of the 
unreflected Brownian motion X. However, for u <0 this cannot be done straight- 
forwardly due to a problem of integrability. In any case, the formal answer given by the 
Girsanov transformation was used to guess the candidate function f of (8) which was 
used for the rigorous proof above. | 


Remark 2. On setting $ = 0 in (7), we see that L(t) is exponentially distributed with 
mean (1 — e~ °)/2u ifu # 0 or Cifu = 0. The exponential distribution of L(t) can also 
be seen as follows (this argument was told to the author by Jim Pitman). For any s 2 0, if 
L(t)>s, then for t, =inf{t20:L,>s}, L(t)=s+L(t°@,) where 0. is the shift 
operator on paths: 0,Z = Z(t + -). Since Z(t,) = 0, it then follows by the strong Markov 
property of Z that 


P(L(t)>t+s|L(t)>s)=P(L(t)>t)  forallt 20. 
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Thus, L(t) is an exponential random variable. To verify the mean of L(t), observe that 
E([Z(tat)) =wE[t at] + E[L(tat)], 
and so on letting t — œ we obtain 


E[L(t)] = E[Z(t)] — wE[t] = C — wE[t]. 


Thus, when u = 0, we obtain E[L(t)] = C, and when u # 0, we can deduce the value of 
E[L(t)] provided we have the value of E[t] as given in Theorem 2 below. One can check 
that the formula for E[L(t)] obtained in this way agrees with that given above. 
Alternative derivations of the first line of (7) can be obtained from [7], Exercises 
(VI.2.10) and (VI.4.12). 

By the strong Markov property of Z, t’ is independent of (L(t), t) and so (5) may be 
rewritten: 


p- (zurea -ELOLE EELDE 
[t] + E[t’] 
(13) 
2 , \2 2 
GEL + AEE + ELI WELD?) / arp, 
(E[t] + Elt)? 


By differentiating (7) suitably and letting a, 8 —> 0, we can find all moments in the above 
involving t and L(t). Note in particular that the existence of finite limits of the 
differentiated expressions as «œ, —0 implies existence of the associated moments 
(cf. Feller [4], p. 435). Since the algebra of these calculations is quite long, Mathematica 
(Version 1.2) on a SUN 3/60 was used to perform it. The supplier of this software 
package is Wolfram Research and a comprehensive reference book is [8]. Now, t’ is 
equal in law to the t for a reflected Brownian motion with drift — u. Thus, the first and 
second moments of t’ can be obtained by replacing u by — u in the corresponding 
expressions for t. The formulas thus obtained for the constituent moments appearing in 
(13) are displayed in the following theorem. 


Theorem 2. Foru=0, 
E(t] = E[t] = C? 
4 


SC 
E[t*] = E[(t’)*] = 3 


E[L(t)] = C 


E{(L(t))*] = 2C? 


3 


5C 
E[{tL(t)] = ~ 
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For u + 0, 
— 2uC) — 1 + 24C 
prs POAT 
2u? 
Fi ee ere 
2u* 
exp(2uC) — 1 — 2uC 
a 
u 
E[(t’)?] = exp(4uC) + exp(2uC) — 6uC exp(2uC) + 2u’?C? — 2 
2u* 
1 — exp( — 2uC) 
Baoe aa 
2u 
(1 — exp( — 2uC))? 
E{(L(t))] = Se 
u 
— 2uC) — C — 2uC) — _4 
aioe a a A ss ced Cn cE ol 


2u? 

When the above expressions are substituted into the right member of (13), after 
simplification, one obtains the expressions in the right member of (6). Thus, we obtain 
the same result (6) as Berger and Whitt ([2], (4.14)). 


Remark 1. The above expressions for E[t] and E[t?] when u = — 1, agree with 
those that can be obtained from the expressions for E[Toọc], Var(Toc) given in Corollary 
3.4.1 of Abate and Whitt [1] where Toc = T. 


Remark 2. Note that if one replaces 0 by — 0 in (6) and one simplifies by multiply- 
ing the resulting numerator and denominator by exp(20C), then one again has formula 
(6). Hence, I’ is the same for u as for — u. This confirms the fact that one gets the same 
formula for I if L(t) in (5) is replaced by U(t + T^) (cf. (4.14) of [2]). 


Acknowledgements 


The author would like to thank Ward Whitt for bringing this problem to her attention 
and for helpful comments on a preliminary draft of this paper. Jim Pitman is also 
thanked for his simple argument that L(t) 1s exponentially distributed. Finally, the 
author thanks the referee for several helpful comments. 


References 


[1] ABATE, J. AND WHITT, W. (1988) Transient behavior of the M/M/1 queue via Laplace 
transforms. Adv. Appl. Prob. 20, 145-178. 

[2] BERGER, A. W. AND WHITT, W. (1992) The Brownian approximation for rate-control throttles 
and the G/G/1/C queue. Dynamic Discrete Event Systems: Theory and Applications 2, 7-60. 


This content downloaded from 
63.135.161.6 on Sun, 01 Sep 2024 03:29:41 UTC 
All use subject to https://about.jstor.org/terms 


1002 R. J. WILLIAMS 


[3] CHUNG, K. L. AND WILLIAMS, R. J. (1990) Introduction to Stochastic Integrations, 2nd edn. 
Birkhauser, Boston. 

[4] FELLER, W. (1971) An Introduction to Probability Theory and Its Applications, Volume 2. Wiley, 
New York. 

[5] Harrison, J. M. (1985) Brownian Motion and Stochastic Flow Systems. Wiley, New York. 

[6] Lions, P. L. AND SZNITMAN, A. S. (1984) Stochastic differential equations with reflecting 
boundary conditions. Comm. Pure Appl. Math. 37, 511-537. 

[7] Revuz, D. AND Yor, M. (1991) Continuous Martingales and Brownian Motion. Springer- 
Verlag, New York. 

[8] WOLFRAM, S. (1991) Mathematica — A System for Doing Mathematics by Computer, 2nd edn. 
Addison-Wesley, Reading, MA. 


This content downloaded from 
63.135.161.6 on Sun, 01 Sep 2024 03:29:41 UTC 
All use subject to https://about.jstor.org/terms 


